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Introduction: 


The present work includes some of the author's original researches on the integer 
solutions of equations and linear systems: 

1. The notion of “general integer solution" of a linear equation with two unknowns 
is extended to linear equations with n unknowns and then, to linear systems. 

2. The properties of the general integer solution are determined (both of a linear 
equation and of a linear system). 

3. Seven original integer algorithms (two for linear equations and five for linear 
systems) are presented. The algorithms are carefully demonstrated and an 
example for each of them is presented. These algorithms can be easily introduced 
into computer. 


INTEGER SOLUTIONS OF LINEAR EQUATIONS 


Definitions and properties of the integer solutions of linear equations. 
Consider the following linear equation: 
(1) Sa x; =b, 

with all a, z 0 à b in Z. 


Again, let h EN, and f, Zr >Z, i= Ln. (Ln means: all integers from 1 to n). 


Definition 1. 


x =x", i=1,n, isa particular integer solution of equation (1), if all x? €Z and 


n 
0 
a =b. 
i-l 





Definition 2. 
x, = f (ky... k,), i=1,n, is the general integer solution of equation (1) if: 
a) Y af (kk) mb; (kk) eZ", 
i=1 
b) For any particular integer solution of equation (1),x, = x? ,i= 1,2, there exist 
(kl, kp) EZ" such that x? = f(k°,... k?) for all i= Ln {i. e. any particular integer 
solution can be extracted from the general integer solution by parameterization}. 


We will further see that the general integer solution can be expressed by linear 
functions. 


h 
For 1X i En we consider the functions f, = ek, +d, with all c;, d; eZ. 


j=l 
Definition 3. 
A — (cj), ; is the matrix associated with the general solution of equation (1). 


Definition 4. 
The integers k,,...,k,, 1<s<h are independent if all the corresponding column 
vectors of matrix A are linearly independent. 


Definition 5. 


An integer solution is s-times undetermined if the maximal number of 
independent parameters is 5 . 


Theorem 1. The general integer solution of equation (1) is (n — 1)-times 
undetermined. 


Proof: 
We suppose that the particular integer solution is of the form: 


D x=} uP +v, izln,withall u,, v, €Z, 
e=1 


P, are parameters of Z , while a<r<n-1. 
Let (x),..., x?) be a general integer solution of equation (1) (we are not interested in 
the case when the equation does not have an integer solution). The solution: 





des j=l,n-1 


a T 
is undetermined (n — ia (it can be easily checked that the order of the associated 


matrix is n —1). Hence, there are n—1 undetermined solutions. Let's consider, in the 
general case, a solution be undetermined (n — 1)-times: 


x = 9 ck; +d,, i-ln with all Cj d; EZ. 


jj 


Consider the case when b = 0. 
Then 


Y'a, -0 

i-l 
It follows: 

= Sa Sect z Ya Zok j o3 =0 
= i=l 
For k, 20, j=1,n—1 it follows that Sad, = 
i=1 

For k, =1 and k, =0, j # ją, it follows that Sac, =0. 


i=1 
Let’s consider the homogenous linear system of n equations with n unknowns: 


Nee. j2lLhn-1 
i=l 


Y ud, =0 
i=1 


which, obviously, has the solution x,=a,, i= Ln different from the trivial one. Hence 








the determinant of the system is zero, i.e., the vectors c; - (c psa j2Lln- 


D- (dodi ) are linearly dependent. 


But the solution being (n — 1)-times undetermined it shows that c pl Ln-1are 


linearly independent. Then (Costo) determines a free sub-module Z of order n—1 in 
Z, of solutions for the given equation. 
Let’s see what can we obtain from (2). We have: 


n n r 
0= Fax = Ya [Sur + v) : 
j=l i=1 e=l 


As above, we obtain: 
n r 
day, = 0 and Jau, =0 
i=l e=l 


similarly, the vectors U, = (tinsa) are linearly independent, h = Lr,U no h= Lr are 
V = (v, V, ) particular integer solutions of the homogenous linear equation. 

Sub-case (a1) 

U,h=1,r are linearly dependent. This gives {U,,...,U,}= the free sub-module of 
order r in Z” of solutions of the equation. Hence, there are solutions from {V,,..., V, } 
which are not from (U,,...,U,] ; this contradicts the fact that (2) is the general integer 


solution. 
Sub-case (a2) 
U, h- Lr, V are linearly independent. Then {U Geel +V is a linear variety 





of the dimension < n—-1 = dim(V,,..., V, ,] and the conclusion can be similarly drawn. 


Consider the case when b + 0. So, Jak =b. 
i=l 


Then: 


n n-l n-l n n 
da 5T + 4j - Sac, + Lad =b; V(k,,...k,,)€Z"". 
i- J= i= 


j=l S isl 





As in the previous case, we obtain Yad, = b and dics -0, Vjzln-1. 
i-l i-l 
The vectors c j^ CE ji: j =1,n-1, are linearly independent because the solution is 
undetermined (n — 1)-times. 
Conversely, if c,,...,c, ,, D (where D = (dyed, J) were linearly dependent, it 

n-l 

would mean that D= ds jC; with all s j scalar; it would also mean that 
j=l 


d dad j da (Ese i Ss (Zae) =0. 
t= i=l j=l 1 


j= i=l 
This is impossible. 


(3) Then {c,,...,¢,,}+D is a linear variety. 


Let us see what we can obtain from (2). We have: 


n n r r n n 
b= SOX = Ya, (Eur + J = S| Sau.) IL Say, 
i=l i=l e=l 1\i=1 i=l 


and, similarly: day, =b and S ai; -0,  Ve=lr, respectively. The vectors 
i=l i-l 
U, = (is-su), e= Lr are linearly independent because the solution is undetermined 
r -times. 
A procedure like that applied in (3) shows that U,..,U,, V are linearly 


independent, where V =(v,,...,.v,). Then {U,,....U,}+V = a linear variety = free sub- 
module of order r « n—1. That is, we can find vectors from [c,,..,c, ,] - D which are 
not from {U ng a +V , contradicting the “general” characteristic of the integer number 


solution. Hence, the general integer solution is undetermined (n - 1)-times. 


Theorem 2. The general integer solution of the homogeneous linear equation 


San: = 0 (all a, ¢Z\{0} ) can be written under the form: 
i-l 


n-l = 
4° week, i=in 
j=l 


(with d, =...=d, =0). 

Definition 6. This is called the standard form of the general integer solution of a 
homogeneous linear equation. 

Proof: 

We consider the general integer solution under the form: 


n-l 
x = NB td, i-ln 
j=l 


with not all d, 20. We'll show that it can be written under the form (4). The 
homogeneous equation has the trivial solution x,=0, i= Ln. There is 
n-l RR 
(p... p, jez” such that $c, p? +d, 20, Vi=l,n. 
j=l 
Substituting: P, =k, +p, j=1,n—1 in the form shown at the beginning of the 


demonstration, we will obtain form (4). We have to mention that the substitution does not 





diminish the degree of generality as P, eZ o k, €Z because j=1,n-1. 


Theorem 3. The general integer solution of a non-homogeneous linear equation is 
equal to the general integer solution of its associated homogeneous linear equation plus 
any particular integer solution of the non-homogeneous linear equation. 

Proof: 


Let's consider that x, = E i-ln,is the general integer solution of the 


j=l 


ij “j? 


associated homogeneous linear equation and, again, let x, = v, i- 1n, be a particular 
n-l 


integer solution of the non-homogeneous linear equation. Then x, = co ¿Wp i= Ln, 
j=l 


is the general integer solution of the non-homogeneous linear equation. 


Actually, Sa -Ya Sica.) Ya (Su Ix =b 
i=1 


i-l 





* 0 à- " 2 
if x,=x,, i-ln, is a particular integer solution of the non-homogeneous linear 
equation, then x, 2 x; -v,, i-Ln, is a particular integer solution of the homogeneous 


linear equation: hence, there is (kj,..., k? ,) eZ" ! such that 
Yok = x° -w, Vi= Ln ; 


i.e.: 


which was to be proven. 


Theorem 4. If x, =>% ck, i=l,n is the general integer solution of a 
jj g 2 
j=1 


homogeneous linear equation CRE ~1 Vj=l,n-1. 
The demonstration is done by reduction ad absurdum. If 3j, 1< j, <Sn—1 such 


that (c; , cs )7 d, *tl,then c, —c,d, with [cross s Vi-ln. 


Cui 


But x; =; , i-ln, represents a particular integer solution as 
0 


n n n 
2 89 = 206. = 1/d; i 2. dé. i 0 
i=l i-l i=] 


(because x, = c; 





p i-Ln is a particular integer solution from the general integer 


solution by introducing k, =1 and k,=0, jz jo. But the particular integer solution 


x, = îs i -ln, cannot be obtained by introducing integer number parameters (as it 


should) from the general integer solution, as from the linear system of n equations and 
n—1 unknowns, which is compatible. We obtain: 
x, = 2 cy key d; k, c i=1,n. 


Y J ijo “Jo Jo Ej 


Leaving aside the last equation — which is a linear combination of other n—1 
equations — a Kramerian system is obtained, as follows: 








Therefore the assumption is false (end of demonstration). 


Theorem 5. Considering the equation (1) with (q,,...,a,)~1, b=0 and the 
n-1 

general integer solution x, = = ck 
j=l 


Ro i = |n, then 
(dif oasis do (Cres cas Vizln. 


Proof: 
The demonstration is done by double divisibility. 


Let's consider ij, 1€ i, <n arbitrary but fixed. x, "En jk; - Consider the 

equation Sax, —-a,x, . We have shown that x,=c,, i= ei is a particular integer 
itig 

solution irrespective of j, a<j<n-1. 


The equation Don = —&, c, ; obviously, has the integer solution x; 2 c;, iz ij. 
izio 


Then (às oro.) divides —a,c,,; as we have assumed, it follows that 


(a,,...,a,) - 1, and it follows that Lares as ab aid), irrespective of j. Hence 


(yaoi ee Oastei) Vi-ln, and the divisibility in one sense was 
proven. 

Inverse divisibility: 

Let us suppose the contrary and consider that di el,n for which 
(A E | od nde dea m we have considered d,, and d,, 
without restricting the generality. d;,ld,, according to the first part of the 
demonstration. Hence, 3d €Z such that d, , =d-d,,, |d |z1. 





il? 


-1 


Sep Cd j -d: diy XA j > 


j=l 


icy E -a,x Zan oce d eh 


i-l izi izi 


= 


where (pes) -]. 


The non-homogeneous linear equation S dde = —a, d; , has the integer solution 


izi 


because a,d,, is divisible by (ao: Let's consider that x, =x , izi, 


is its particular integer solution. It follows that the equation Xar = 0 the particular 
i-l 


solution x, 2 x^, iz i, x, = d, , which is written as (5). We'll show that (5) cannot be 


obtained from the general solution by integer number parameters: 


n-1 

Dek =x, izi, 
(6) ee n-1 

drd, Ok di 


But the equation (6) does not have an integer solution because d-d,,|d,, thus, 
contradicting, the “general” characteristic of the integer solution. 
As a conclusion we can write: 


Theorem 6. Let’s consider the homogeneous linear equation Sax, =0, with all 
i=1 


a, €Z\{0} and (a,...,a,) - 1. 
h — 
Let x, = Ok i-ln, with all c; eZ, all k, integer parameters and let's 
j=l 


consider h e N be a general integer solution of the equation. Then, 
1) the solution is undetermined (n — 1)-times; 


2) Vj-Ln-1l we have (cose) sale 


3) Vi=Ln we have (c,,..0, 1) ~ (apa 0,4, ). 


il?***? "in-l 


The proof results from theorems 1,4 and 5. 


Note 1. The only equation of the form (1) that is undetermined n -times is the 
trivial equation 0-x, t... 0: x, — 0. 


Note 2. The converse of theorem 6 is not true. 


Counterexample: 
x, --k-4 k, 
(7) x, = 5k, +3k, 
x,=7k,—k,; k,k eZ 
is not the general integer solution of the equation 
(8) —13x, + 3x, — 4x, =0 


although the solution (7) verifies the points 1), 2) and 3) of theorem 6. (1, 7, 2) is the 
particular integer solution of (8) but cannot be obtained by introducing integer number 
parameters in (7) because from 


-k + k, -1 
5k, +3k, =7 
Jk =k =2 


it follows that k = g Z and k= i € Z (unique roots). 


REFERENCE 
[1] Smarandache, Florentin — Whole number solution of linear equations and 


systems — diploma thesis work, 1979, University of Craiova (under the 
supervision of Assoc. Prof. Dr. Alexandru Dincá) 
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AN INTEGER NUMBER ALGORITHM TO SOLVE LINEAR 
EQUATIONS 


An algorithm is given that ascertains whether a linear equation has integer number 
solutions or not; if it does, the general integer solution is determined. 


Input 
A linear equation a,x,+...+a,x, =b, with a,beZ, x, being integer number 


unknowns, i — Ln „and not all a; =0. 


Output 
Decision on the integer solution of this equation; and if the equation has solutions 
in Z, its general solution is obtained. 


Method 
Step 1. Calculate d —(a,,..., a, ). 


Step 2. If d /b then “the equation has integer solution"; go on to Step 3. If d Xb 
then “the equation does not have integer solution”; stop. 
Step 3. Consider h=1. If |d|z1, divide the equation by d; consider 


a, -a;/d, i2 ln, bi bld. 


Step 4. Calculate a — min 








a,| and determine an i such that a; =a . 


Step 5. If a + 1 then go to Step 7. 
Step 6. If a=1, then: 
(A) x =-lax Pau Pen X UL Poe tb), 
(B) Substitute the value of x, in the values of the other determined 


unknowns. 
(C) Substitute integer number parameters for all the variables of the 
unknown values in the right term: k,,k,,...,k,.,, and k, , 


respectively. 
(D) Write, for your records, the general solution thus determined; stop. 
Step7. Write down all a,, j zi and under the form: 


a= ng tT 


a; b 
b=aq+r where q, = Sas „al 


Step 8. Write x, = qux, —..—q, X; í,7 Qi 7 q,X, + 9 — by. Substitute the 
value of x, in the values of the other determined unknowns. 
Step 9. Consider 


11 


a; = =d; 
Gu Shu b=r 
i f and 
Ciu > Tin X; > t, 
h=h+1 
a, =, 


and go back to Step 4. 


Lemma 1. The previous algorithm is finite. 
Proof: 
Let's a,x, +...+a,x, =b be the initial linear equation, with not all a; = 0; check 


for min|a, | =a,#1 (if not, it is renumbered). Following the algorithm, once we pass 
a, #0 


for 





from this initial equation to a new equation: a,x, ^ ax, ...-- a,x, =b', with lai < la; 





i=2,n, |b'|<|b| and a, 2-a,. 


It follows that min|a,| < min |a, |. We continue similarly and after a finite number 
azo! * azi | * 


of steps we obtain, at Step 4, a:=1 (the actual a is always smaller than the previous a, 
according to the previous note) and in this case the algorithm terminates. 


Lemma 2. Let the linear equation be: 
(25) ax,+a,x,+..+4,x, =b, with min |a, | =a, and the equation 
az 





(26) -at thx% +.. +r xX, =r, with 1t,-—x,-4,x,—-..—q4,x,*q, where 
. ^4. š a; b 0 
r-a-àq;, i-2n, r=b-aq while q;=|—|, r=|—|. Then x=x, 
a a, 


X,—X)..,X,-— X» is a particular solution of equation (25) if and only if 





t =t =x -qX —...— q, X? +q, X.X, =X? is a particular solution of equation 
(26). 


0 
n? 


x 2X, XX, =x, is a particular solution of equation (25) 


a,x? + a,x} + a, x? -be ax? + G, + aq) Fart (ri + aq, )x. -aqctre 





0 0 0 0 0 E 0 0 0 
nX, t.e +r X, Q(X, = pă e ON, +a) HTS al rA rX Sr 


n n 


e t 21), x, 2 x), x, 2 x5 is a particular solution of equation (26). 


k,,+4d,, i= 1n, is the general solution of equation 


l 


Lemma 3. x, =c, k *...*c 
(25) if and only if 
(28) t — 0G bye +-+ a, Rei Cna Copa EFC. JR, = 
—(d, +q,d, +...+4,d,)+a, 


in-l 


X= Caka tC sk, td; jz2n 


jn-l 
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is a general solution for equation (26). 
Proof: 


0 0 0 0 0 0 ` E z 
t =t =X X 0 aX, td, X = X5,..,X, =X, is a particular solution of 





the equation (25) o x, =x), x, = 32,.., x, =x? is a particular solution of equation (26) 
€» dk, =k) €Z,...,k, =k) €Z such that 

x, =c k? +..4¢, k? +d, =x, i=ln e3k =k eZ, k, =k? eZ, 
such that 

X UK +...+C 


in-l 
0 ow PE 
ia ei +d, Tz X; , l= 2,n > 
and 
0 0 
t = (Cy, + q,0, * nen dK, ied (er F GaG E dC), — 





(d, +q,d, +...+q,d,)+ q =-X} — 49x; -—q,X +q =t; 


Lemma 4. The linear equation 
(29) ax+a,x, +...+a,x, =b with la,| = 1 has the general solution: 
x =—-(a,k, +...+.a,k, — b)a, 
(30) x,=k, eZ 
i=2,n 
Proof: 
Let's consider x, 2 xj, x, 2 x2,..., x, = x? , a particular solution of equation (29). 
JE =x}, k, =x? , such that x, = (—a,25 +a -b)a, IUE ENT EN 


Lemma 5. Let's consider the linear equation ax, +a,x, t... a,x, =b, with 








min|d, —a, and a,=a,g,, i=2,n. 
a, 20 


a Ss 


Then, the general solution of the equation is: 
x} = -(q,k, Teck q,k, = q) 
x -keZ 
i0 
Proof: 
Dividing the equation by a, the conditions of Lemma 4 are met. 


Theorem of Correctness. The preceding algorithm calculates correctly the 
general solution of the linear equation a, x, t... a,x, =b, with not all a; 2 0. 

Proof: 

The algorithm is finite according to Lemma 1. The correctness of steps 1, 2, and 3 
is obvious. At step 4 there is always min |a, | as not all a; 2 0. The correctness of sub- 


step 6 A) results from Lemmas 4 and 5, respectively. This algorithm represents a method 
of obtaining the general solution of the initial equation by means of the general solutions 
of the linear equation obtained after the algorithm was followed several times (according 


13 


to Lemmas 2 and 3); from Lemma 3, it follows that to obtain the general solution of the 
initial linear equation is equivalent to calculate the general solution of an equation at step 
6 A), equation whose general solution is given in algorithm (according to Lemmas 4 and 
5). The Theorem of correctness has been fully proven. 


Note. At step 4 of the algorithm we consider a := min la,| such that the number of 
a, 


iterations is as small as possible. The algorithm works if we consider a = la, # max 
s=l,n 








a Ss 


but it takes longer. The algorithm can be introduced into a computer program. 


Application 
Calculate the integer solution of the equation: 
6x, —12x, —8x, + 22x, =14. 


Solution 

The previous algorithm is applied. 

1. (6,—12,—8,22)= 2 

2. 2114 therefore the solution of the equation is in Z. 
3. h :- 5; [2| z 1; dividing the equation by 2 we obtain: 
3x, = 6x, — 4x, *11x, 2 7. 


























4. a= min[B|, |-6|, |-4], [11] 23. i=1 
5. asl 
7. -6z3.(-2)*0 
-4z3.(-2)*2 
11=3-3+2 
723.241 
8. x =2x, +2x, 3X1 2-1, 
9. 
a, =0 a, =—3 
a, =2 b=1 
a,=2 X, m 
h=2 


4. We have a new equation: 
-3t —0-x, +2x, +2x,=1 














a= min {|-3 325 2) and 
LD-3 
5. az1 
7. -3=2-(-2)4+1 
0=2-0+0 
2=2-1+0 
1=2-0+0 


14 


8. x, =2t,+0-x, — x, -0— t, . Substituting the value of x, in the value 
determined for x, we obtain: x, = 2x, — 5x, +31, - 2t, +2 


a, =0 X, =t; 
h=3 
4. We have obtained the equation: 
1-t,+0-x, —2-t,+0-x,=1, a=1, and i=1 
6. (A) t, 2 -(0: x, -2t, +0-x, - D-12265 +1 
(B) Substituting the value of 7, in the values of x, and x, previously 
determined, we obtain: 
x, =2x, —5x,+4t, +5 and 
beta Oe +3t, +2 
(C) x, =k, x,=k,, tp =k, k,k,,k, eZ 
(D) The general solution of the initial equation is: 
x, = 2k, - 5k, - Ak, +5 


x, =k, 
x, = —k, + 3k, +2 
x, =k, 


k,,k,,k, are parameters eZ 


REFERENCE 

[1] Smarandache, Florentin — Whole number solution of equations and 
systems of equations — part of the diploma thesis, University of Craiova, 
1979. 
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ANOTHER INTEGER NUMBER ALGORITHM TO SOLVE 
LINEAR EQUATIONS (USING CONGRUENCES) 


In this section is presented a new integer number algorithm for linear equation. 
This algorithm is more “rapid” than W. Sierpinski's presented in [1] in the sense that it 
reaches the general solution after a smaller number of iterations. Its correctness will be 
thoroughly demonstrated. 


ANOTHER INTEGER NUMBER ALGORITHM TO SOLVE LINEAR 
EQUATIONS 


Let's us consider the equation (1); (the case a,b €Q, i=1,n is reduced to the 
case (1) by reducing to the same denominator and eliminating the denominators). Let 
d = (a,,...,a,). If dib then the equation does not have integer solutions, while if d |b 
the equation has integer solutions (according to a well-known theorem from the number 
theory). 

If the equation has solutions and dz we divide the equation byd. Then d=1 
(we do not make any restriction if we consider the maximal co-divisor positive). 

Also, 

(a) If all a, the equation is trivial; it has the general integer solution 

x, =k, €Z, i=1,n, when 5-0 (the only case when the general solution is 


n -times undetermined) and does not have solution whenb + 0. 
(b) If 3i, 1<i<n such that a,=+1 then the general integer solution is: 


Pu \ 
ee MU a and x, =k, €Z, se{l,....n}\{i} 
j=l 


jzi 
The proof of this assertion was given in [4]. All these cases are trivial, therefore 
we will leave them aside. The following algorithm can be written: 


Input 
A linear equation: 


Q) Janbao bez a *tli-ln, 
i=l 


with not all a; =0 and (a,,...„a,)=1. 


Output 

The integer general solution of the equation. 
Method 

l.h=l,p=l 
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2. Calculate min {|r 


I<i, j<n 








, r=a,(moda,), 





r|< la il} and determine r and the pair 


(i, j) for which this minimum can be obtained (when there are more possibilities we have 
to choose one of them). 
3. If |r| # go to step 4. 


If |r| 2 1, then 








(A) Substitute the values thus determined of these unknowns in all the 
statements (p), p = 1,2,... (if possible). 


(B) From the last relation (p) obtained in the algorithm substitute in all 
relations: (p — 1), (p — 2 ys(D 

(C) Every statement, starting in order from (p B 1) should be applied the 
same procedure as in (B): then (p — 2) respectively. 





(D) Write the values of the unknowns x, i=1,n, from the initial 


equation (writing the corresponding integer number parameters from 
the right term of these unknowns with K,...,k, ,), STOP. 





4. Write the statement (p): x; =t,- AT 
a. 
J 
5. Assign Xj, h=h+1 
a, =r p=ptl 


The other coefficients and variables remain unchanged go back to step 2. 
The Correctness of the Algorithm 


Let us consider linear equation (2). Under these conditions, the following 
properties exist: 


Lemma 1. The set M =Í |r 


Proof: 

Obviously M c N' and M is finite because the equation has a finite number of 
coefficients: n , and considering all the possible combinations of these, by twos, there is 
the maximum AR? (arranged with repetition) = n^ elements. 





, r=a,(moda,), 0<|r|<|a,|} has a minimum. 
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Let us show, by reduction ad absurdum, that M + Ø. 
M#@ = a; 2 0(moda;) Vi, j = Ln. Hence a, = O(moda,) Vi, j= 1,n. Or this 


Ps Ln, which is equivalent to 





is possible only when =|a, 








a. 


(a,,...4,)=4a,, Vi eln. But (a,.,a,)-1 are a restriction from the assumption. It 
follows that la, = Ln, Vi el,n a fact which contradicts the other restrictions of the 


assumption. 
M z 0 and finite, it follows that M has a minimum. 


a;l, viel,n. 








Lemma 2. If |r| = min M , then Ir|« 
<i,jsn 
Proof: 


We assume conversely, that Ji), 1 € ij € n such that |r] > a, | ; 


Then r|> min aj] - la; | 1. ISj,Sn.Leta,, 1< py <n, such that la, |» a; and 


1<j<n 


a, is not divided by a. 


There is a coefficient in the equation, 





a; | which is the minimum and the coefficients 
are not equal among themselves (conversely, it would mean that (q,,..,a,)=a,=+1 
which is against the hypothesis and, again, of the coefficients whose absolute value is 
greater that la;, | not all can be divided by a, (conversely, it would similarly result in 
(a,,..,a,) = d; * +1. 


We write |a, /a,, ] =q, EZ (integer portion), and r=a, —q a, €Z. We have 


" „| € |r|. Thus, we have found an 7, which |r;| « |r| 





a, = rn(moda, ) and 0 «|| «Ja, | <a, 
contradicts the definition of minimum given to |r]. 


Thus|r|<|a,|, vi eLn. 








a; 


Lemma 3. If |r| = min M = 1 for the pair of indices (i, j), then: 
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is the general integer solution of equation (2). 
Proof: 
Let x,=x°?, e- Ln, be a particular integer solution of equation (2). Then 
ai 





=f 
x; €Z (because a; — r = Ma;) such 


dk, =x? eZ, se Iss fi, j} and f, =x} + 
a 


s) 


that: 











and x, =k, 2x,, se(L.., n] MA j}. 


Lemma 4. Let [r| and (i, j) be the pair of indices for which this minimum can 
be obtained. Again, let's consider the system of linear equations: 


n 
at, t TX, + 3 a,x, =b 


(3) seli.j) 





Then x, 2x2, e= Ln isa particular integer solution for (2) if and only if x, = x°, 
a-r 





e €(1,...,n) Mj] and t, =t, =x; + x, is the particular integer solution of (3). 
j 


Proof: 
x, =x, e=1,n is a particular solution for (2) if and only if 


n n 
a, -r 
bax be» ps sata [t l e)re 





a. 


e=l s=l j 





n 
a,—r 
0 0 0 0 0 0 0 
€» af, +x, + > a,x, =b and t, =x,+——-x, eZ Sx =X 


s-l d; 


sei, j} 


ee f, Ss n) \ uj and f, =f, is a particular integer solution for (3). 


Lemma 5. The previous algorithm is finite. 
Proof: 
When |r|=1 the algorithm stops at step 3. We will discuss the case when 


r| EN’. We will show that the row of 





|r|z 1. According to the definition of r, 
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r—s successively obtained by following the algorithm several times is 
decreasing with cycle, and each cycle is not equal to the previous, by 1. Let 7, be 


the first obtained by following the algorithm one time. Ul #1 then go to step 4, 
, Viz ln. 

Now we shall follow the algorithm a second time, but this time for an 
equation in which 7, (according to step 5) is substituted by a,. Again, according 











and then step 5. According to lemma 2, |r| « |a, 


to lemma 2, the new |r| written |r,| will have the propriety: |z| « ||. We will get 
to |r|= 1 because |r| 2 1 and |r| « oc, and if |r| z 1, following the algorithm once 
again we get |r| « Inl and so on. Hence, the algorithm has a finite number of 
repetitions. 


Theorem of Correctness. The previous algorithm calculates the general 
solution of the linear equation correctly (2). 

Proof: 

According to lemma 5 the algorithm is finite. From lemma 1 it follows 
that the set M has a minimum, hence step 2 of the algorithm has meaning. When 
|r| 2 1 it was shown in lemma 3 that step 3 of the algorithm calculates the general 


integer solution of the respective equation correctly the equation that appears at 
step 3). In lemma 4 it is shown that if |r| z 1the substitutions steps 4 and 5 


introduced in the initial equation, the general integer solution remains unchanged. 
That is, we pass from the initial equation to a linear system having the same 
general solution as the initial equation. The variable A is a counter of the newly 
introduced variables, which are used to successively decompose the system in 
systems of two linear equations. The variable pis a counter of the substitutions of 
variables (the relations, at a given moment between certain variables). 

When the initial equation was decomposed to |r| = 1, we had to proceed in 
the reverse way, i.e. to compose its general integer solution. This reverse way is 
directed by the sub-steps 3(A), 3(B) and 3(C). The sub-step 3(D) has only an 
aesthetic role, i.e., to have the general solution under the form: x, = f,(k,,...,k,_,), 


i-ln, f, being linear functions with integer number of coefficients. This “if 


possible" shows that substitutions are not always possible. But when they are we 
must make all possible substitutions. 


Note 1. The previous algorithm can be easily introduced into a computer 
program. 


Note 2. The previous algorithm is more “rapid” than that of W. 
Sierpinski's [1], i.e., the general integer solution is reached after a smaller number 
of iterations (or, at least, the same) for any linear equation (2). 

In the first place, both methods aim at obtaining the coefficient +1 for at least 
one unknown variable. While Sierpinski started only by chance, decomposing the 
greatest coefficient in the module (writing it under the form of a sum between a 
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multiple of the following smaller coefficient (in the module) and the rest), in our 
algorithm this decomposition is not accidental but always seeks the smallest |r| 
and also choose the coefficients a, and a; for which this minimum is achieved. 
That is, we test from the beginning the shortest way to the general integer 
solution. Sierpinski does not attempt to find the shortest way; he knows that his 
method will take him to the general integer solution of the equation and is not 
interested in how long it will take. However, when an algorithm is introduced into 
a computer program it is preferable that the process time should be as short as 
possible. 


Example 1. 

Let us solve in Z^ the equation 17x — 7y + 10z = -12. 
We apply the former algorithm. 

1. h=1,p=1 

2.r=3,i=3,j=2 

3. [3| 21 go on to step 4. 





10 — 
4. (1) y=4,- u-htz 
5. Assign 
y=t, h:=2 
4,73 p=2 


with the other coefficients and variables remaining unchanged, go back to 
step 2. 
2. r=-]l, i=l, j=3 
AES 
x = —l(-3t, — (771) - 12) = 3t, — 74, — 2 
17-(- -1- 
(-1) T 1—17 
3 


(A) We substitute the values of x and z thus determined into the only 
statement (p) we have: 


(1 yHt4+z=—-—17t, +43t — 72 








s (m« 71): ( 12) =17t, + 421, — 72 


(B) The substitution is not possible. 
(C) The substitution is not possible. 
(D) The general integer solution of the equation is: 


x = 3k, — Tk, +12 
y =-17k, + 43k, — 72 
z=—-17k,+42k,-72; k,k eZ 
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[1] 
[2] 
[3] 
[4] 
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INTEGER NUMBER SOLUTIONS OF LINEAR SYSTEMS 


Definitions and Properties of the Integer Solution of a Linear System 


Let's consider 


(y Saab isa 
j=l 





a linear system with all coefficients being integer numbers (the case with rational 
coefficients is reduced to the same). 





Definition 1. x, =X j=l,n, is a particular integer solution of (1) if 





y J 


n 

0 . 0 E 

x; eZ, j=1,n and > a,x, =b, izlm. 
j=l 


Let's consider the functions f,:Z" > Z, j= In, where h eN'. 


Definition 2. x pH fi (ks kr) J = Ln , is the general integer solution for (1) if: 


(a) DATE ess) -b, i=1,m, irrespective of (k,,...k, eZ; 


jel 

(b) Irrespective of x, = x, „j=Ina particular integer solution of (1) there is 
(kf, kp) e Zsuch that. f, (kr, kp) =x), j= Ln. (In other words the general solution 
that comprises all the other solutions.) 


Property 1. 

A general solution of a linear system of m equations with n unknowns, 
r(A) =m «n, is undetermined (n —m)-times. 

Proof: 

We assume by reduction ad absurdum that it is of order r, 1€ r En — m (the case 
r = 0,i.e., when the solution is particular, is trivial). It follows that the general solution is 
of the form: 


X Hp, cup, +y 
(Si) 
X, =U Pi fe FU P. +V, Up VieZ 
p, = parameters e Z 
We prove that the solution is undetermined (n — m) -times. 


The homogeneous linear system (1), resolved in r has the solution: 
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x, = ty tet xx, 
p k 
Xp = 5 PME d Pon 
Let x, = x°, i- Ln, be a particular solution of the linear system (1). 
Considering 
Xm+1 = D ` ku 
x,=D-k, 
we obtain the solution 
uS 1 0 
x = Dpa Kna teet D, se +x 
m m 0 
Xn = D, ` ka Tec D, "k, tX, 
eta 0 
Xma > D: ka + Xml 
x, 2 D-k, +x, k, = parameters e Z 


which depends on the n — m independent parameters, for the system (1). Let the solution 
be undetermined (n —m)-times: 


X = Caki Pal AE a td 


n-m 


S 
(S2) kit. WEN AOL 


nn—-m n-m 


x, =C 


nl 
Cj, d, € Z, k, = parameters € Z 

(There are such solutions, we have proved it before.) Let the system be: 
aX taux, = b 


At FF co aux, = D, 
b eZ. 


I. The case b, 0, i=1,m results in a homogenous linear system: 


x, =unknowns eZ, d;, 
aX; t ..ba,x,-0;i- Lm. 

(Sj — —aj (ck * e, uk us td) t au (cuf + + eu uk, us td, )= 0 
0 = (acu Toc a, Lai Ted (6G). Se Oe a + (a.d, Tec a, d, ) 
Vk; EZ 

For k =...=k,_,, =0=> 4a,d,+...+a,d, =0. 
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For, 2. ky = kpa ==, = 0 and k,=1 > 


n-m 


> (a;c;, Tod aC, + (a.d, Tuck a, di" )- 0 





QC, +...+4,,C,, 20, Vil m,Vhzln-m. 


in "nh 


The vectors 


fei 





Can 
are the particular solutions of the system. 
V,, h=1,n-m also linearly independent because the solution is undetermined 
(n-m)-times {V,,...,V,_,,}+d is a linear variety that includes the solutions of the system 


obtained from (Sz). 
Similarly for (S1) we deduce that 


(U;,) 





are particular solutions of the given system and are linearly independent, because (S1) is 


(V. 
undetermined (n —m)-times, and V 2| | isa solution of the given system. 


V, 


Case (a) U,,...,U,, V= linearly dependent, it follows that {U,,...,U,}is a free 
sub-module of order r «n- m of solutions of the given system, then, it follows that there 
are solutions that belong to{V,,...,V,_,,}+d and which do not belong to {U,,...,U,, a 


fact which contradicts the assumption that (S1) is the general solution. 


Case (b) U,,...,.U, , V = linearly independent. 
{U 1525 U. ; +V isa linear variety that comprises the solutions of the given system, which 
were obtained from (S1). It follows that the solution belongs to {V,,...,V,_,,} +d and does 
not belong to {U eee Ul 2 +V ,a fact which is a contradiction to the assumption that (S1) is 
the general solution. 

II. When there is an i €1,m with b, + 0 then non-homogeneous linear system 

A,X, +..+ apă, -b,i- lm 

(So) > aj (ek +--+ Cp ku + ) tas (ouk ete 
it follows that 

= (ac, Jack AC za (ainm Fase ua + (a.d, fub a, d, )= b, 


in“ nn-m 


+d, )- b, 


nn—m kam 
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For kj =...=k,_,, =0 >a,d,+..+a,d, =b; 
For k 2..-k,,-k =k, =0 and k, =1> 





fel, > 


— (anc ECOLE )* (aud, t... a; d, ) =b, it follows that 


j in "nj in^n 


ae -0 E 





in “nj -— 


2 
a,d,+...+a,,d, =b; 
(ea) 


V,=|: |,j=ln-m, are linearly independent because the solution (S2) is 


Či 


undetermined (n — m )-times. 
(d) 
Q) — V, jeln-m,and E | 
d, 
are linearly independent. 
We assume that they are not linearly independent. It follows that 
(Sic +--+ S, Lucus 
d= sV +t Sp ,V, 47 
Lp, a ah ae So 
Irrespective of i= l,m 
b = aud, + --+ ad, = du (sic, Teck 


in^n 


n-mCIn-m )+ et ain ics tec SEL er Je 


= (ac, Tec GC Si Ft (CAT s teat Gia. )o =0. 
Then, b, =0, irrespective of i= Lm , contradicts the hypothesis (that there is an i e Lm, 
b, #0). It follows that V,,...,V,_,,,d are linearly independent. 

asas) +d is a linear variety that contains the solutions of the non- 
homogeneous system, solutions obtained from (S2). Similarly it follows that 
{G.G} +V is a linear variety containing the solutions of the non-homogeneous 


system, obtained from (S). 
n -m > r it follows that there are solutions of the system that belong to 





“21!” means “to prove that" 


(V... V, ,j* d and which do not belong to {G,,...,G,}+V , this contradicts the fact that 


(S1) is the general solution. Then, it shows that the general solution depends on the n — m 
independent parameters. 


Theorem 1. The general solution of a non-homogeneous linear system is equal to 
the general solution of an associated linear system plus a particular solution of the non- 
homogeneous system. 

Proof: 
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Let's consider the homogeneous linear solution: 
aX cua x =0 


, (AX 20) 
Q, X Fe FAX, = 0 
with the general solution: 
X = Cu ++ Cin kn td 


X, = Cy Ky FF Can mn S T d, 


and 


E 
XQ = A 


0 
X, =X, 
with the general solution a particular solution of the non-homogeneous linear system 
AX =b; 
= 0 
X = Ok, Hach ES A + 
(2!) 
x, 2 Cy ky toe, Kn + dy X) 


nn-m n—m 


is a solution of the non-homogeneous linear system. 


We note: 
Fawn Was 


SML 


(vector of dimension m ), 


AX = A(Ck+d+x° )=A(Ck+d)+ AX" =b+0=b 
We will prove that irrespective of 


_ 0 
x, = y 


x, = Ya 
there is a particular solution of the non-homogeneous system 
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Kpm = ade € Z 
with the property: 


m 0 0 0. .0 
x, = Cu Het Cn knn rd t =), 


0 0 0 0 
k, coc, mKn-m tdt m y, 


nn—-m n-m 


x, =C 


(0) 


We note Y? i | 
0 
Yn 


We'll prove that those k? €Z, j=1,n-m are those for which A(cx gs d)- 0 


(there are such X? €Z because 


nl 


x,20 


x, =0 
is a particular solution of the homogeneous linear system and X = CK + d is a general 
solution of the non-homogeneous linear system) 
A(CK" d « X* -Y*) - A(CK* +d)+ AX" AY =0+b-b=0 


Property 2 The general solution of the homogeneous linear system can be written 
under the form: 





(SG) 
Xi = e dave Le 
(2) 
X, = Caki ape LE 
k, is a parameter that belongs to Z (with d, =d, =...=d,=0). 
Poof: 


(SG) = general solution. It results that (SG) is undetermined (n — m) -times. 
Let's consider that (SG) is of the form 

X= Cy Py t.t Cin D, +H 
(3) 

X, > CAD, tC +d, 
with not all d; =0; we'll prove that it can be written under the form (2); the system has 
the trivial solution 


m-mPn-m 


28 


x, =0eZ 
x,-0eZ 
it results that there are p, €Z, j = Ln-m, 
X =P H.+ CnaP n td =0 
(4) 
X, = CaP? H.+ CnP na +d, =0 
Substituting p, 2 k, + p^, j=1,n-m in (3) 
k, eZ 
p eZ 


J 


lon eZ, 


peZ 0 

p =z —k,-p,-p,eZ 

which means that that they do not make any restrictions. 
It results that 


0 0 
x = cuki H.+ Co en + Cupi tec Cin m Dum + d,) 


X, cO POET (cap? +46), Pop + d,) 
But 
CaP? uo CV VP, td, =0, h= Ln (from (4)). 
Then the general solution is of the form: 
X = Cuk +-+ Cin-nK, 


n-m 


X, = C,íK, +--+ 0, SK, 


nn-m n-m 





k, = parameters eZ, j = l,n- m; it results that d, = d, =..=d,=0. 


Theorem 2. Let's consider the homogeneous linear system: 
A,X, +...+4,,x, =0 
, 


QUT taux, = 0 


r(A)-m, CANE T) E 1, h=1,m and the general solution 
X, = Cki +--+ e, VK 


n—m 


X, = Chi bs Ck 


nn-m n-m 


then 
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(615 Gica Gao im) 
irrespective of h = 1,m and i=1,n. 
Proof: 


Let's consider some arbitrary h €1,m and some arbitrary i €1,n; 
aX +-+ UL Xr dax 








on Tru UL, oU, X 


i? 


Because 





[d 69d a aues) a, 


it results that 
d = (dijs ii Ahiri Ary ) 


irrespective of the value of x, in the vector of particular solutions. 


Xi 








For k, = k, = ...= k, ,, = 0 and k, «1 we obtain the particular solution: 


Xp = Cu 


x,=c, — dle, 


X, = C, 
For k,=k,=...=k,_,,;=0 and k,,-1 it results the following particular 
solution: 
X, m Cin n 
Xi = Ciim > d | Cim > 
X = Cun-m 
hence 
dlc, j-ln-m- | CON RR S 
Theorem 3. 
If 


X = Cak +... + 6, Kk 


n-m 


X, = Caki FF C, mK 


k, = parameters €Z, c; eZ being given, is the general solution of the homogeneous 
linear system 


aX +.. +a, X, = 0 
, r(A)=m<n 
CX Tob Apă, = 0 
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then REN =1, Vi=ln-m. 
Proof: 
We assume, by reduction ad absurdum, that there is j, €l,n-m: (4, deut ) =d 


5 njo 





we consider the maximal co-divisor > 0; we reduce to the case when the maximal co- 
divisor is —d to the case when it is equal to d (non restrictive hypothesis); then the 
general solution can be written under the form: 


X, = Cuki Feet Cup dk, Ee uk 


n-m 
(5) 
Xp = Caki +-+ Cpp dk, OL VEL. 
where d aaan; Cy = Cy and (eae ad 
We prove that 
x, = Cy, 
X, = Ci, 


is a particular solution of the homogeneous linear system. 


We'll note: 
(k, ) 
f c, T Cn. ds cus) 
Cai : , k= 


nl] '** Cj, ase Cnm 


k; 
k, =R 


x=C-k the general solution. 


(agan) 
We know that AX = 0 > A(CK)=0, A=|: 


a a 


nl*** “mn 

We assume that the principal variables are x,,...,x,, (if not, we have to renumber). 
It follows that x, ,,...., x, are the secondary variables. 
For k,=..=k, ,-k,,-..-k,,-0 and k, =1 we obtain a particular 
(c, d) fe) (c,, \ x, = Cy; 


—0-2A!: —-d.A|: —A: =0> 





solution of the system 
Tp 


X, = €, d c, d ome c 


njo njo n njo 
is the particular solution of the system. 
We'll prove that this particular solution cannot be obtained by 


3l 


Ky = yk Feet Cup dk, c.c, uk = rj, 
(6) 

X, = Caki t.c dk, t. cu uk um Co, 

Xn = Coah Tec C, dk; Tec NA RI = m+l jy 
(7) 


X, = Cy hy bee dk, Feo o, LK 


nn-m “n-m  . um 


Cul $ Cuaj ut CoL a-m 
C itd XC 
hj nj n,n—m 
SS kes =—¢€¢Z 
Jo d 
Ca Sia Cui, exe Col n-m 
Cha Cnj d C, n-m 


(because d #1). 
It is important to point out the fact that those k, = $ , j=1,n— m , that satisfy the 


system (7) also satisfy the system (6), because, otherwise (6) would not satisfy the 
definition of the solution of a linear system of equations (i.e., considering the system (7) 
the hypothesis was not restrictive). From X, €Z follows that (6) is not the general 


solution of the homogeneous linear system contrary to the hypothesis); then 





(e = 1, irrespective of j = l,n— m. 


Property 3. Let's consider the linear system 
aX, o aux, =, 


X ba yx = bn 
b, eZ, r(A)=m<n, x, = unknowns EZ 
Resolved in R , we obtain 


x 2 fi (Xs) 


d; 


: , Xj- X, are the main variables, 
Xm = Ín (Cate) 
where f, are linear functions of the form: 
+.. + cix, +e, 
d, 


t 


i 
f. = Co mal 
= 





> 





where ci 


mj? 


d,e eZ; i-lm,j-ln-m. 


32 


e; ; ; NET ; : 
If — e Z irrespective of i=1,m then the linear system has integer solution. 
Proof: 
For 1X i &m, x, eZ, then f, eZ. Let's consider 


X m = Usi 
X, = u,k, 
e 

EN 1 €, 

X = Vo Ka Het V, KE, 
1 

Dom k rS mp i Em 

Xn — Vol m+l t V, n 


m 


a solution, where u,,, is the maximal co-divisor of the denominators of the fractions 
i 
mj 


d, 


t 


C 








,i-lm,j-ln-m calculated after their complete simplification. 


i 
i __ m+j m+j 


m+j — d 


t 


v €Z is a (n— m)-times undetermined solution which depends on 


" k,) but is not a general solution. 


mar? 


n — m independent parameters (k 


Property 4. Under the conditions of property 3, if there is an 


Wu ; ' e, y, E. e; 
i elm:f,-upjXx,4T-ctux,*—- with up; eZ, j=ln—m, and — eZ then the 
io io 
system does not have integer solution. 
Proof: 


Vx X, in Z, it results that x, ¢Z. 


ml?** 
Theorem 4. Let's consider the linear system 
aX +..+4,,%, =D, 


Qs +-+ Gy x, =D, 


ab; EZ, x;> unknowns €Z, r(A)=m<n. If there are indices 1<i, <...<i,, € n, 


ij? 


i, €{1,2,...n}, h= 1,m, with the property: 
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A= #0 and 
mi, mi, 
b, a, a 
A, = is divided by A 
4 
b, a ni, I a ni 
dy = Au, b, 
Ax, = : 1 | is divided by A 
a mi pia A mina m 


then the system has integer number solutions. 
Proof: 
We use property 3 

d, =A, i=1,m; e, = Ay, ,h=l1,m 


Note 1. It is not true in the reverse case. 


Consequence 1. Any homogeneous linear system has integer number solutions 
(besides the trivial one); r(A) 2 m «n. 
Proof: 
A, =0:A, irrespective of h = Lm. 


in 


Consequence 2. If A = +1, it follows that the linear system has integer number 
solutions. 
Proof: 


Ax : (£L), irrespective of h=1,m; 
A, EZ. 


ln 
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FIVE INTEGER NUMBER ALGORITHMS 
TO SOLVE LINEAR SYSTEMS 


This section further extends the results obtained in chapters 4 and 5 (from linear 
equation to linear systems). Each algorithm is thoroughly proved and then an example is 
given. 

Five integer number algorithms to solve linear systems are further given. 


Algorithm 1. (Method of Substitution) 

(Although simple, this algorithm requires complex computations but is, 
nevertheless, easy to implement into a computer program). 

Some integer number equation are initially solved (which is usually simpler) by 
means of one of the algorithms 4 or 5. (If there is an equation of the system which does 
not have integer systems, then the integer system does not have integer systems, then 
Stop.) The general integer solution of the equation will depend on n — 1 integer number 
parameters (see [5]): 


= pO GO ai jaiai 
(p,) X; =f (k^. sy > b= ln; 
where all functions i are linear and have integer number coefficients. 


This general integer number system (p,) is introduced into the other m-1 
equations of the system. We obtain a new system of m —1 equations with n — 1 unknown 
variables: 

ks i, =1,n-1, 
which is also to be solved with integer numbers. The procedure is similar. Solving a new 
equation, we obtain its general integer solution: 

(p) IOS (Ok) ee nat, 
where all functions y are linear, with integer number coefficients. (If, along this 
algorithm we come across an equation, which does not have integer solutions, then the 
initial system does not have integer solution. Stop.) 

In the case that all solved equations had integer system at step (j), ls j<r 
( r being of the same rank as the matrix associated to the system) then: 

G-D 2 eO (LU) DĂ = 
(p) ki (c i: i,2ln-j-l, 
fi are linear functions and have integer number coefficients. 

Finally, after r steps, and if all solved equations had integer solutions, we obtain 
the integer solution of one equation with n —r+ 1 unknown variables. 

The system will have integer solutions if and only if in this last equation has 
integer solutions. 

If it does, let its general integer solution be: 

(p) KEP EUR S i =l,n-r+1, 
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where all f” are linear functions with integer number coefficients. 
We'll present now the reverse procedure as follows. 
We introduce the values of ki”, i, =1,n—r+1, at step p, in the values of 


ke „i =l,n-r+2 
from step (p, ,). 

It follows: 

ich OP ce (RO) =e Rae), 


(r-1) 


,=Ln-r-—1, from which it follows that g; are linear functions with integer 


uber coefficients. 
Then follows those (p, ;) in (p, ,) and so on, until we introduce the values 
obtained at step (p,) in those from the step (p,). 


It will follow: 
= gl (kt... Xe) 


notation g, (k, xs i =l,n, with all g, most obviously, linear functions with 


integer number coefficients (the notation was made for simplicity and aesthetical 
aspects). This is, thus, the general integer solution, of the initial system. 


The correctness of Algorithm 1. 

The algorithm is finite because it has r steps on the forward way and r — 1 steps 
on the reverse, (r « +00). Obviously, if one equation of one system does not have (integer 
number) solutions then the system does not have solutions either. 

Writing Sfor the initial system and S, the system resulted from step (p;), 


1< j<r-2, it follows that passing from (p,) to (p,,,) we pass from a system S, to a 


system S,,, equivalent from the point of view of the integer number solution, i.e. 
Ci 2 ln TI 
k^ -f,d-Lln-jtl, 
which is a particular integer solution of the system S, if and only if 


Kk 


ESI 


is a particular Du solution of the system S,,, where 


1) 
=a. (as AEQ ij, = Ln- J. 


Hence, their 2 integer solutions are also equivalent (considering these 
substitutions). Such that, in the end, resolving the initial system S is equivalent with 
solving the equation (of the system consisting of one equation) S, , with integer number 
coefficients. It follows that the system S has integer number solution if and only if the 
systems S, have integer number solution, 1< j S r- 1. 


E , plug J, 


jH 





Example 1. By means of algorithm 1, let us calculate the integer number solution 
of the following system: 
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5x—-T7y-2zt6w-6 
(S) | 


—4x + 6y- 3z+1lw=0 


Solution: We solve the first integer number equation. We obtain the general 
solution (see [4] or [5]): 

x-t +24, 
y-h 

(p) 
z= -t,+5t, + 3t, -3 
w=t, 

where t, t, t, EZ. 
Substituting in the second, we obtain the system: 


(S,) St, — 23t, +24 +9=0. 

Solving this integer equation we obtain its general integer solution: 
t, — K 

(P2) t, =k, +2k,+1 


t, =9k,+23k,+7 

where k,, k, eZ. 

The reverse way. Substituting (p,) in (p,) we obtain: 

X = 3k, +4k, +2 
y-k 
z= 31k, + 79k, +23 
w — 9k, + 23k, +7 

where k,, k, € Z, which is the general integer solution of the initial system (S). Stop. 


Algorithm 2. 
Input 
A linear system (1) without all a; 20. 


Output 
We decide on the possibility of an integer solution of this system. If it is possible, 
we obtain its general integer solution. 


Method 
1.t=1,h=1,p=1 
2. (A) Divide each equation by the largest co-divisor of the coefficients of the 
unknown variables. If you do not obtain an integer quotient for at least one equation, then 
the system does not have integer solutions. Stop. 
(B) If there is an inequality in the system, then the system does not have integer 
solutions. Stop. 
(C) If repetition of more equations occurs, keep one and if an equation is an 
identity, remove it from the system. 
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3. If there is (i, j,) such that 





d, i | —] then obtain the value of the variable 
x, from equation i; statement (T,). 

Substitute this statement (where possible) in the other equations of the system and 
in the statement (T, ,), (H,) and (P,) for all i, ^, and p. Consider 1: 1&1, remove 
equation (i,) from the system. If there is no such a pair, go to step 5. 


4. Does the system (left) have at least one unknown variable? If it does, consider 
the new data and go on to step 2. If it does not, write the general integer solution of the 
system substituting k,, k,,... for all variables from the right term of each expression 


which gives the value of the unknowns of the initial system. Stop. 
5. Calculate 


and determine the indices i, ji, ją as well as the r for which this minimum can be 
calculated. (If there are more variables, choose one arbitrarily.) 


a = min [|rja, = r(moda, ), 0<|r|< 
ijoj Ji J2 








i, 


l a; —r 
6. Write: x, =t, ———x 
J2 : 
ij 


possible in all the equations of the system and in the statements (T,), (H,) and (P,) for 


statement (H,). Substitute this statement (where 


i? 


all t, h, and p. 
7. (A)If az, consider x, =1,, h=h+1, and go on to step 2. 
(B) If a = 1, then obtain the value of x, from the equation (i); statement 
(P,). Substitute this statement (where possible in the other equations of the 
system and in the relations (7,), (H,) and (P, ,) for all t, ^, and p. 


Remove the equation (i) from the system. 
Consider À:- h 1, p := p+1, and go back to step 4. 


The correctness of algorithm 2. Let consider system (1). 


Lemma 1. We consider the algorithm at step 5. Also, let 
M -[r TN INE 
Then M z Ó. 
Proof: 
Obviously, M is finite and M c N . Then, M has a minimum if and only if 
M + 0. We suppose, conversely, that M = 0. Then 
di, = O(moda,, Vi jej 





; d = r(moda, ), 0<|r|< 








ay 


It follows as well that 
a; = O(moda, ) VAs Siders 
That is 








la, |- OV joja 
We consider an i, arbitrary but fixed. It is clear that 


ay, 
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(apioa n) | 4; #0, Vj 


(because the algorithm has passed through the sub-steps 2(B) and 2(C). But, because it 
has also passed through step 3, it follows that 








a, | * L Vj, 
but as it previously passed through step 2(A), it would result that 
a, | =b Vj. 








This contradiction shows that the assumption is false. 


Lemma 2. Let's consider a, , = r(moda, ) . Substitute 





dij, 
in system (A) obtaining system (B). Then 
x, = 2 j=l,n 
is the particular integer solution of (A) if and only if 
x, =x; j * j, and t, =x; — 


is the particular integer solution of (B). 


Lemma 3. Let a, + and a, obtained at step 5. 
Then 0 «a, <a, 
Proof: 


It is sufficient to show that a, <|a,|, V i, j because in order to get a, , step 6 is 








d 
obligatory, when the coefficients if the new system are calculated, a, being equal to a 
coefficient form the new system (equality of modules), the coefficient on (i, j,). 


Ea. 





Let a, ,, with the property la, ü 


i Let a,, with 











Hence, a, 2 = min (ja, d. ; | > ; there is such an element 











a. j Hi 


,j=l,n 








because la; " | is the minimum of the coefficients in the module and not all a, ; 


are equal (conversely, it would result that (a, ,,....,a,,) ~ a; ;, Vj elr, the algorithm 


passed through sub-step 2(A) has simplified each equation by the maximal co-divisor of 


its coefficients; hence, it would follow that =], Vj= Ln , Which, again, cannot be 








d; j 
real because the algorithm also passed through step 3). Out of the coefficients a,, we 


choose one with the property a,, # Ma, there is such an element (contrary, it would 





result (a, ,...„4,,)-la,. | but the algorithm has also passed through step 2(A) and it 
19J tn 10Jm p 


would mean that la; n | —] which contradicts step 3 through which the algorithm has also 


passed). 
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Considering qə = la, is, / dli, NI eZ and r-a,, -qqa,, EZ, we have 


10J59 


= n(moda,, ) and 0< Ins 





< <a. We have, thus, obtained an 7, with 








d, in | S la, i 


lo I So 








In] < a, which is in contradiction with the very definition of a,. Thus a, « a;|, Vi, j. 

Lemma 4. Algorithm 2 is finite. 

Proof: 

The functioning of the algorithm is meant to transform a linear system of m 
equations and n unknowns into one of m, x n, with m, <m, n, <n, thus, successively 
into a final linear equation with n—r+1 unknowns (where r is the rank of the 
associated matrix). This equation is solved by means of the same algorithm (which works 
as [5]). The general integer solution of the system will depend on the n—1 integer 
number independent parameters (see [6] — similar properties can be established also the 
general integer solution of the linear system). The reduction of equations occurs at steps 
2, 3 and sub-step 7(B). Step 2 and 3 are obvious and, hence, trivial; they can reduce the 
equation of the system (or even put an end to it) but only under particular conditions. The 
most important case finds its solution at step 7(B), which always reduces one equation of 
the system. As the number of equations is finite we come to solve a single integer number 
equation. We also have to show that the transfer from one system m, xn, to another 


m; Xn,,, is made in a finite interval of time: by steps 5 and 6 permanent substitution of 


variables are made until we to a=1(we to a=1 because, according to lemma 3, all 
a — s are positive integer numbers and form a strictly decreasing row). 


Theorem of correctness. 

Algorithm 2 correctly calculates the general integer solution of the linear system. 

Proof: 

Algorithm 2 is finite according to lemma 4. Steps 2 and 3 are obvious (see also 
[4], [5]). Their part is to simplify the calculations as much as possible. Step 4 tests the 
finality of the algorithm; the substitution with the parameters k, k,,... has 
systematization and aesthetic reasons. The variables 7, A, p are counter variables (started 
at step 1) and they are meant to count the statement of the type T, H,P (numbering 
required by the substitutions at steps 3, 6 and sub-step 7(B); h also counts the new 
(auxiliary) variables introduced in the moment of decomposition of the system. The 
substitution from step 6 does not affect the general integer solution of the system (it 
follows from lemma 2). Lemma 1 shows that at step 5 there is always a, because 
ØM cN“. 

The algorithm performs the transformation of a system m, xn, into another 
m xn equivalent to it, preserving the general solution (taking into account, 
however, the substitutions) (see also lemma 2). 


i+1? 


Example 2. Calculate the integer solution of: 
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12x - Ty + 9z =12 
— 5y+ 8z+10w=0 
0z+ Ow=0 

15x t21z-469w = 3 


Solution: 
We apply algorithm 2 (we purposely selected an example to be passed through all 
the steps of this algorithm): 
1. t=1,h=1,p=1 
2. (A) The fourth equation becomes 5x + 7z+23w=1 
(B)- 
(C) Equation 3 is removed. 
3. No; go on to step5. 
.a=2 and i=1, j 22, j, - 3, and r=2. 
6. z- t, +y, the statement (H,). Substituting it in the 
12x-2y * 9t, E 
3y+9t,+10w =0 
5x +7y+7t,+23w =1 


CA 


. a#1 consider z = t, := 2 , and go back to step 2. 


. No. Step 5. 
.a-landiz2,j,-4,j,-2,and r=1. 





6. y « t, — 3w, the statement (H,). Substituting in the system: 
—12x * 2t, +9t -6w -12 
3t, +8t + w=0 
5x +7t,+7t +2w=1 
Substituting it in statement (H,), we obtain: 
z=1 +t, — 3w, statement (H,). 
7. w=-3t, — 8t, statement (P). 
Substituting it in the system, we obtain: 
—12x — 20t, + 574 212 
| 5x +t,- 9t =1 
Substituting it in the other statements, we obtain: 
z=10t, +25t, (H,)" 
y=10t,+24t,, (H,)" 
h:=3, p:=2 , and go back to step 4. 
4. Yes. 
2.— 
3. t, = 1—- 5x *9t,, statement (T;). 
Substituting it (where possible) we obtain: 
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{-1 12x+237t, = —8 (the new system); 
z= 10-50x+1154,, (Hj)" 
y= 10-50x+114t,, (H;)" 
w=-34+15x+ 35, (B)' 
Consider f = 2 go on to step 4. 
4. Yes. Go back to step 2. (From now on algorithm 2 works similarly with that 
from [5], with the only difference that the substitution must also be made in the 
statements obtained up to this point). 
2.— 
3. No. Go on to step 5. 
5. a= 13 (one three) and i=1, j = 2, j, 2 1, and r «13. 
6. x « t, - 2t, , statement (H,). 
After substituting we obtain: 
41120, + 131, 2 -8 (the system) 
z=10-50f +154, (H,)"; 
y =10— 504 +14t, (H,)"; 
w=-3+15t, -5t, (2) 
t, = 1- 5t, — f (m); 
7. x 5 t4, h= 4 and go on to step 2. 
2.— 
3. No, go on to step 5. 
5. a=5 and i=1, j =1,j,=2 andr=5 
6. t, =t, +9t,, statement (H,). 
Substituting it, we obtain : 
5t, + 13t, = —8 (the system). 
z= 104851, +15, (Hj); 
y= 10-761, -14t,, (H,)" ; 
x= 19t + 2t, (H,)'; 
w--3-30t- 5t, (B)'ó 
t= l-M&- 4 (ID 
7. t — th 5 and go back to step 2. 
2.— 
3. No. Step 5. 
5. a=2 and i=1, j =2, j, =1 and r=-2. 
6. t, — t, — 3t, statement (H,). After substituting, we obtain: 
5t; — 2t, = —8 (the system). 
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z=104+85t,—240t,, (H); 
y=10+76t,—214t,, (H,)’; 
y= 19t, -555,. (EL); 
w=-3-30t,+85t,, (B); 
t --1-1454415,, T)"; 


t= 9t,+26t,,  (H,)5 
7. t, = t,,h = 6 and go back to step 2. 
2, 
3. No. Step 5. 


5. a=1 and i=1, j =2, j,,r=1. 
6. t, =t,+2t, statement (H,). After substituting, we obtain: 
t, —2t, =—8 (the system) 
z 210-3951. —2404, (H,)"; 
y =10— 3924, — 2144, (H,)"; 
x= —91t,-55t,, (H); 
w=-3+1401,+851,,  (P)'; 
t= 1+68t, +41, (7); 
t=  -—435-264, (H)* 
t, = -5t 3ti (Huy 
7. t, = 21, — 8 statement (P,). Substituting it in the system we obtain: 0=0. 


Substituting it in the other statements, it follows: 
z = —10304 + 3170 


y = —918t, + 2826 
x =—237t, + 128 
w= 3651, — 1123 
t, = 1774, —543 

t, 21121, +344 
134, + 40 
Sig 16 


statements of no importance. 


f 
f 


I 


Consider h = 7, p = 3, and go back to step 4. t €Z 
4. No. The general integer solution of the system is: 
x-—2XIK, + 128 


y = —918k, + 2826 
z=1030k, +3170 
w= 365k, -1123 
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where k, is an integer number parameter. 
Stop. 


Algorithm 3. 

Input 

A linear system (1) 

Output 

We decide on the possibility of an integer solution of this system. If it is possible, 
we obtain its general integer solution. 


Method 




















1. Solve the system in R”. If it does not have solutions in IR", it does not have 
solutions in Z” either. Stop. 

2 f=1,t=1,h=1,g=1 

3. Write the value of each main variable x, under the form: 


(E a) xm 2 4,5, +q, [xw = /A, 
J J 


<|A,|, A; #0, |r| «|A, 


the right term are integer number variables: either of the secondary variables of the 
system or other new variables introduced with the algorithm). For all i, we write 
Lo = Ke. 


Uli 











with all q;, q;, 7, r, A;in Z such that all r, 

















lj 





(where all x, of 


4. (E ie ;) 33/77 -nY, tn-0 where (v A ;) are auxiliary integer number 
J 


variables. We remove all the equations (F ; ; ) which are identities. 
5. Does at least one equation (F, ;) exist? If it does not, write the general 


integer solution of the system substituting k,,k,,... for all the variables from the right 


term of each expression representing the value of the initial unknowns of the system. 
Stop. 


6. (A) Divide each equation (F, ;) by the maximal co-divisor of the coefficients 
of their unknowns. If the quotient is not an integer number for at least one i, the system 
does not have integer solutions. Stop. 

(B) Simplify —as in m - all the fractions from the statements (E A iJ: 


7. Does r,, exist having the absolute value 1? If it does not, go on to step 8. If it 


Jo 
does, find the value of x, from the equation (F, i y write (7,) for this statement, and 
substitute it (where it is possible) in the statements (£ * ) (7) (H ^ ) (c, ) for all 
i, t, h and g. Remove the equation (F, ^ ) Consider f:= f *1, t: t1, and go back 


to step3. 
8. Calculate 
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I} 


hj 





a = min {|r 
ijoj 





„7, =r(modr, ), 0< Ir« 


and determine the indices i,, j, j, as well as the r for which this minimum can be 


obtained. (When there are more variables, choose only one). 


; d cR. ; 
9. (A) Write x, —z,— M Xx „> Where z, is a new integer variable; statement 


(H,). 


(B) Substitute the letter (where possible) in the statements 


(£,.) Ea) ©). Q1,.). (G, ) forall i, t, h and g. 


(C) Consider h:=h+1. 
10. (A) If a+ go back to step 4. 


(B) If a2 1 calculate the value of the variable x, from the equation (F n J 
relation (c; ) Substitute it (where possible) in the statements (E | ) (7, ) (H h ) (c...) 





for all i, t, h, and g. Remove the equation (F, ) Consider g := g +1, f :=f+1 and go 
back to step 3. 


The correctness of algorithm 3 





j*; A, (with all Fj, To Ap n, My being 
jen 


Lemma 5. Let i be fixed. Then p Fx r) 


integers, n, < n,, A, £0 and all xX; being integer variables) can have integer values if 


and only if (5. xod A) eee 
Proof: 
The fraction from the lemma can have integer values if and only if there is a 


z EZ such that 
Ny n; 
Dr tr IA o ze Pnx,-Aztr-0, 
j=m j=m 


which is a linear equation. This equation has integer solution <> (r, saci A,) |r. 


Lemma 6. The algorithm is finite. It is true. The algorithm can stop at steps 1,5 or 
sub-steps 6(A). (It rarely stops at step 1). 

One equation after another are gradually eliminated at step 4 and especially 7 and 
10(B) (F, | ;) - the number of equation is finite. 


If at steps 4 and 7 the elimination of equations may occur only in special cases 
elimination of equations at sub step 10 (B) is always true because, through steps 8 and 
9 we get to a=1 (see [5]) or even lemma 4 (from the correctness of algorithm 2). 
Hence, the algorithm is finite. 


Theorem of Correctness. 
The algorithm 3 correctly calculates the general integer solution of the system (1). 
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Proof: 

The algorithm if finite according to lemma 6. It is obvious that the system does 
not have solution in R” it does not have in Z” either, because Z” c R” (step 1). 

The variables f, 7, h, g are counter variables and are meant to number the 


























statements of the type E, F, T, H and G , respectively. They are used to distinguish 
between the statements and make the necessary substitutions (step 2). 

Step 3 is obvious. All coefficients of the unknowns being integers, each main 
variable x; will be written: 


which can assume the form and conditions required in this step. 
Step 4 is obtained from 3 by writing each fraction equal to an integer variable 
Y, ; (this being x, €Z). 


Step 5 is very close to the end. If there is no fraction among all (E A ;) it means 


that all main variables x, already have values in Z , while the secondary variables of 
the system can be arbitrary in Z, or can be obtained from the statements T, H or G 
(but these have only integer expressions because of their definition and because only 
integer substitutions are made). The second assertion of this step is meant to 
systematize the parameters and renumber; it could be left out but aesthetic reasons 
dictate its presence. According to lemma 5 the step 6(A) is correct. (If a fraction 
depending on certain parameters (integer variables) cannot have values in Z , then the 
main variable which has in the value of its expression such a fraction cannot have 
values in Z either; hence, the system does not have integer system). This 6(A) also 
has a simplifying role. The correctness of step 7, trivial as it is, also results from [4] 
and the steps 8-10 from [5] or even from algorithm 2 (lemma 4). 


Ther initial system is equivalent to the “system” from step 3 (in fact, (E f. ;) as 


well, can be considered a system) Therefore, the general integer solution is preserved 
(the changes of variables do not prejudice it (see [4], [5], and also lemma 2 from the 
correctness of algorithm 2)). From a system m, xn, we form an equivalent system 


with m, «m, and n 


m,, 1% nia i+] i+] 


2. 


<n,. This algorithm works similarly to algorithm 


Example 3. Employing algorithm 3, find an integer solution of the following 
system: 


SX 4x, +22x,- 8x, =25 
6x, + 46x, —12x, =2 
4x, +3x,—x,+9x,; = 26 
Solution 
1. Common resolving in R° it follows: 
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230 


X 
-3 
x,*2x, +24 
a 4 
11,2 
3 











3 X, = pp = (E, ;) 
n= -4n E) 


2x; +3y,-1=0 (F) 
4. x *2x,-4y, =0 (F) 
x,-3y,*2-0 (F) 


5. Yes. 

6. — 

7. Yes: ls —]. Then x, =3y, —2 the statement (7). Substituting it in the 
others, we obtain: 





2x,-1 

x, =-7x,+ 6y,,-4+—4 (E, .) 
x,+6y,,-4 

x, = 6 4 r (E, ») 
3y,, -2+2 

x, = i PA +8 +B (E, .) 


Remove the equation (F, : ) 
Consider f := 2, t := 2; go back to step 3. 





x =-7x, + 6y,,-4+ — : (E,,,) 
3 X = ya 5 + si = (E, z) 
xc 11y,48 (E, .) 
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2x, +3y,-1=0 (F) 


x,*2y,-4y, =0 (Ea) 
5. Yes. 
6. — 
7. Yes |r,| 31. We obtain x, = —2y,, + 4y,, statement (7 ). Substituting it in 
the others we obtain: 


—4 8y, -1 
ME a (E, Ji 


4 Y» + yu *5 (Es) 
x, = -11y +8 (5-4); 


x, =—28y,, + 20y + 





Remove the equation (F, a) 
Consider f := 3, t := 3 and go back to step 3. 


3. 
2y,-2y,-1 
X 7-22), 430, + (Ei) 
X= Yt Vy t5 (E, zi) 
LS —lly,, +8 (£,) 
4. 2y,*2y5 +3y,-1=0 (F) 
5. Yes. 
6. — 
7. No. 
8. a=1 andi, =1, j =31, j, =22,and r=1. 
9. (A) Y» = z Y; (statement (H,). 


(B) Substituting it in the others we obtain: 


Dy. Da dal 
x, 7722 ~ 30g, + 309, 744 Pa (E) 


X7 Xt — Yu +5 (E: 
x, --lly, +8 (5) 


2y +22, + y4 71-0 (Ba) 
X, =—2 y +42 -4y (T,)' 
(C) Consider h := 2 


10. (B) y, = 1- 2y; — 2z, statement (G, ). 
Substituting it in the others we obtain: 
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x =—40y -9273+27 (E) 


x,= 3y,9 3z 44 (E 2) 
x =-l1y; +8 (Ey 
x, = 6y 4122, —4 (T,)" 
Y» =2y,,+3z, -1 (H,)' 


Remove equation (F, |) 
Consider g := 2, f := 4 and go back to step 3. 


3. 
x, = -40y,, - 922, +27 (Eis) 
x,— 3y,+ 3z,+4 (E) 
sity. +8 (Ei) 
4. - 


5. No. The general solution of the initial system is: 
x,=—40k,-92k,+27, from (E,,) 


“= 3k+ 3k, +4, from (E) 
x, -—llk, +8, from (Es) 
x= 6k,+12k, —4, from (T,)" 
x,= 3k, 2r from (7) 


where k,,k, eZ. 


Algorithm 4 

Input 

A linear system (1) with not all a; 2 0. 

Output 

We decide on the possibility of integer solution of this system. If it is possible, we 
obtain its general integer solution. 


Method 

1. hz 1, v=1. 

2. (A) Divide every equation i by the largest co-divisor of the coefficients of the 
unknowns. If the quotient is not an integer for at least one i, then the system 


does not have integer solutions. Stop. 


(B) If there is an inequality in the system, then it does not have integer 
solutions 


(C) In case of repetition, retain only one equation of that kind. 
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(D) Remove all the equations which are identities. 


3. Calculate a= min (la; 
uJ 





a, £ 0) and determine the indices i,, jọ for which this 


minimum can be obtained. (If there are more variables, choose one, at 
random.) 
4. If a#1 goon to step 6. 
If a=1, then: 
(A) Calculate the value of the variable x, from the equation i, note this 


statement (V, ) 


(B) Substitute this statement (where possible) in all the equations of the 
system as well as in the statements (V,_,), (H n for all v and h. 


(C) Remove the equation i, from the system. 


(D) Consider v=v+1. 
5. Does at least one equation exist in the system? 
(A) If it does not, write the general integer solution of the system substituting 
k,, k,,.... for all the variables from the right term of each expression 
representing the value of the initial unknowns of the system. 
(B) If it does, considering the new data, go back to step 2. 
6. Write all a,;, j # jo and b, under the form : 








d = i i,j ti with | < a, | ` 


igj ? 


b j =4, ud, t With |r, 








«la, |. 
7. Write x, = -5 q X; +q, +t,» statement (H,). 
J* Jo 
Substitute (where possible) this statement in all the equations of the system as 
well as in the statement (V.). (H ah for all v and A. 


8. Consider 


Xj m, h:=h+|, 


Aij e Tia J * Jos 


a,. = ta b, = +r, , 


iojo © ~~~ iojo? ~to 0 


and go back to step 2 


The correctness of Algorithm 4 

This algorithm extends the algorithm from [4] (integer solutions of equations to 
integer solutions of linear systems). The algorithm was thoroughly proved in our previous 
article; the present one introduces a new cycle — having as cycling variable the number of 
equations of system — the rest remaining unchanged, hence, the correctness of algorithm 
4 is obvious. 
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Discussion 

1. The counter variables h and v count the statements H and V , respectively, 
differentiating them (to enable the substitutions); 

2. Step 2 ((A)+(B) + (C)) is trivial and is meant to simplify the calculations (as 
algorithm 2); 

3. Sub-step 5 (A) has aesthetic function (as all the algorithms described). 
Everything else has been proved in the previous chapters (see [4], [5], and 
algorithm 2). 


Example 4. Let us use algorithm 4 to calculate the integer solution of the 
following linear system: 


Dai ce gx OX. = —2 
4x,+3x,  -46x,—5x,- 19 
olution 
. hxlvzl 


S 

1 

2: 

3. a=3 andi=1l,j=1 
4. 3+1. Goon to step 6. 
6 


. Then, 
—]23.(-3)42 
6=3-2+0 
-223.0-2 


Te XQ 33 23,1, statement (H ji Substituting it in the second equation we 
obtain: 
At + 3x, +12x, - x, — 5x, =19 
8. x =t, h >= 2, a, = 0, ag = +2, a4,:-0, a, := +3, b=-2. 
Go back to step 2. 
2. The equivalent system was written: 
3t, qx. --2 
4t, +3x, +12x, — x, 5x, = 19 
3. a=], i=2,j=4 
4. 1=1 
(A) Then: x, =4t, +3x, +12x, — 5x, —19 statement (V, ). 
(B) Substituting it in (H Aa we obtain: 
x, =-7t, -6 x,- 21x, +10x, +38, (H,) 
(C) Remove the second equation of the system. 
(D) Consider: v:= 2. 
5. Yes. Go back to step 2. 


2. The equation +34 +2x = —2 is left. 
3. a=2 andi=1,j=3 
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4. 2+2, go to step 6. 
6. 
132 42.2-1 
—2 = +2(-1)+0 
7. x, =—2t, * t, —1 statement (H,). 
Substituting it in (H i ) A (V.). we obtain: 
x, =351,-6 x, —21t, 10x, +59 (H,)" 
x, =—20t, +3x, +12t, -5x,-31 (V,)' 
8. x, =t,,h=3, a):=—-1, a4: 42, b,:=0 , (the others being all = 0). Go back 
to step 2. 
2. The equation —5f, + 2t, = 0 was obtained. 
3. a=l1,and i=1,j=1 
4. 1=1 
(A) Then ¢, = 2t, statement (V,). 
(B) After substitution, we obtain: 
x, =49t, 6x,*10x,*59 (H,)"; 
x, =—28t, +3x, — 5x, 3231 (v)" ; 
x, =—3t, (H,)'; 
(C) Remove the first equation from the system. 
(D) v=3 
5. No. The general integer solution of the initial system is: 
x, = 49k, — 6k, + 10k, + 59 


X, = k, 
x, =—3k, -1 
x, = —28k, + 3k, - 5k, — 31 
X = k, 
where (k,, ,, k,) eZ. 
Stop. 
Algorithm 5 
Input 
A linear system (1) 
Output 


We decide on the possibility of an integer solution of this system. If it is 
possible, we obtain its general integer solution. 
Method 
1. We solve the common system in R”, then it does not have solutions in IR", 
then it does not have solutions in Z” either. Stop. 
2. f=1,v=1, h=1 
3. Write the value of each main variable x, under the form: 
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(E) EG = Dai at [Sos «Ja. 


with all q;, q, rj T, A, from Z such that all |; « 








«|A,|. 





A, # 
(where all d —S of the right term are integer variables: either from the 
secondary variables of the system or the new variables introduced with the 


algorithm). For all i, we write r, = A, 


4. (E, 2 DI T Yi th =O as (,;) are auxiliary integer variables. 
J 


Remove all the equations (F | ;) which are identities. 


5. Does at least one equation (F, ) exist? If it does not, write the general integer 


solution of the system substituting k,, k,,... for all the variables of the right 
number of each expression representing the value of the initial unknowns of 
the system. Stop. 

6. (A) Divide each equation (F, ,) by the largest co-divisor of the coefficients of 


their unknowns. If the quotient is an integer for at least one i, then the system 
does not have integer solutions. Stop. 
(B) Simplify — as previously ((A)) all the functions in the relations (£ jJ: 





7. Calculate a = min |r, 51, # 0} , and determine the indices i), j; for which this 
LJ 


minimum is obtained. 
8. If a#1, go onto step 9. 
If a=1, then: 


(A) Calculate the value of the variable oe from the equation (F, ;) write 
(V.) for this statement. 

(B) Substitute this statement (where possible) in the statement (E f. Ji 
(V, CER (zu, ). for all i, v, and h. 

(C) Remove the equation (E f. 3 

(D) Consider v :2 v +1, f = f +1 and go back to step 3. 


9. Write all r ,, j # jọ and r, under the form: 


Fi 











Hj = A, i di; T n )J ? 


oa. 








4, tn 


hU A; io 


10. (A) Write x, 2 —* q, x; * q, +t, statement (H,). 
J* Jo 


(B) Substitute this statement (where possible) in all the statements (E " i} 


(F, ) C.) (H, ). 
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(C) Consider h:=h+1 and go back to step 4. 


The correctness of the algorithm is obvious. It consists of the first part of 
algorithm 3 and the end part of algorithm 4. Then, steps 1-6 and their correctness were 
discussed in the case of algorithm 3. The situation is similar with steps 7-10. (After 
calculating the real solution in order to calculate the integer solution, we resorted to the 
procedure from 5 and algorithm 5 was obtained). This means that all these insertions 
were proven previously. 

















Example 5 
Using algorithm 5, let us obtain the general integer solution of the system: 
3x, +6x, 42x; =0 
4x, = 2x, —7x, =-—l 
Solution 
1. Solving in IR we obtain: 
m —6x; 2X, 
3 
E —2x, +7x; -1 
4 
2. f=l,v=1, h=1 
—2 
3. (E. ): x, =2x,+ M 
2x, + 3x. -1 
(E2): H= dud EST ci 
4. (F,):-2x,-3y,, =0 
(Fa): 2x, + 3x; — 4y —1=0 
5. Yes 
6. — 
7. i=2 and i, =2, j, =3 
8. 241 
9. 3=2-14+1 
—4 =2-(-2) 
-1=2-0-1 


10. x, =—x,+2y,, +t, statement (H i After substitution: 


23; 





[E EE —2x,—4y,-2t + 


xTM 





(F,)':x 425-120 
Consider h := 2 and go back to step 4. 
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4. (Ei): 2x, — 3 =0 
(Fa) :24 43 —1=0 
. Yes. 


NN 


7. a=1andi,=2, j,=5 
(A) x, =—21, +1 statement (V.) 
(B) Substituting it, we obtain: 


, —2x 
(E,,)": x, =-6t,+2-4y,, 4 — 





(2a): x, =—2t +1+ yy 
(H,)': x, =34+1-1+2y, 
(C) Remove the equation (F2 ) 


(D) Consider v 22, f =2 and go back to step 3. 
-—2X; 


Se ESL X 776, -4y $24 —7* 
(Enz): x; 2 25 + y, +1 

4. (Ei): -2x,-3y, =0 

5. Yes. 

6. — 

7. a=2 and i,=1,j,=4 

8. 2z1 

9, -3--2.(D-1 


10. (A) x, = —y4, +t, statement (H,) 
(B) After substitution, we obtain: 


! = EM 
(P ede = 61, -4y, p 


(eae -yn — 2t, = 0 
Consider h := 3, and go back to step 4. 
4. (E): —yy—2t, =0 
5. Yes 
6. — 
7. a=1 and i, =1, j, =21 (two, one). 
(A) y,, =-2t, statement (V, ). 
(B) After substitution, we obtain: 


(C) Remove the equation (5) ; 
(D) Consider v = 3, f = 3 and go back to step 3. 
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3. (Ez, ): x, 2-65 — 4y — 2t, +2 
(E, ,): x; = 2t», +1 
a 


5. No. The general integer solution of the system is: 


x, =—2k, +k, +1, from (E y; 
x, 23k, +2k, -L from (H,)* 
ae 3k, , from (H,)* 
x; =—2k, +1, from (V, ); 


where (k, k,,k;) eZ. 
Stop. 


Note 1. Algorithm 3, 4, and 5 can be applied in the calculation of the integer 
solution of a linear equation. 


Note 2. The algorithms, because of their form, are easily adapted to a computer 


program. 


Note 3. It is up to the reader to decide on which algorithm to use. Good luck! 
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